NASA 

Technical 

Paper 


2396 


naSA .TP-2396 19850004525 


November 1984 


Algorithm for Astronomical, 
Extended Source, Signal- to- 
Noise Ratio Calculations 

R. R. Jay roe 





NASA 

Technical 

Paper 

2396 

1984 


Algorithm for Astronomical, 
Extended Source, Signal-to- 
Noise Ratio Calculations 


R. R. Jayroe 

George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 


NASA 

National Aeronautics 
and Space Administration 


Scientific and Technical 
Information Branch 





TABLE OF CONTENTS 


Page 


GENERAL DESCRIPTION 1 

USER INSTRUCTIONS 1 

EQUATIONS USED 3 

REFERENCES 7 

APPENDIX A - OUTPUT EXAMPLE 10 

APPENDIX B - PROGRAM LISTING 12 


iii 




TECHNICAL PAPER 


ALGORITHM FOR ASTRONOMICAL, EXTENDED SOURCE, SIGNAL- 
TO-NOISE-RATIO CALCULATIONS 

GENERAL DESCRIPTION 


The software computes signal- to-noise ratios (S/N) as a function of observation 
time and observation times as a function of signal-to-noise ratios for viewing an astro- 
nomical extended source with a telescope and focal plane detector array outside the 
Earth’s atmosphere. The extended source is assumed to be square or rectangular, 
uniform in brightness, and centered on the optic axis of the telescope. A combined 
S/N is computed for all the detector picture elements (pels) affected by the extended 
source, as well as one for a single pel, which may be completely or only partially 
filled (at the user's discretion) by the extended source. The user is also allowed 
to shift the extended source image in any direction by an amount equal to or less 
than the pel dimensions. 

The cosmic background within the field-of-view is also assumed to be uniform 
in brightness. The user may request a monochromatic or polychromatic extended 
source and background. The user also has the choice of requesting a uniform or 
blackbody spectral distribution for the source and background. 

The optical system is characterized by the telescope entrance pupil diameter 
and the system focal ratio. The user can also include optical filters by specifying 
their transmission characteristics. 

The focal plane detector is assumed to be a charge coupled device (CCD) array 
that can have square or rectangular shaped pels and dead spaces between the pels. 
The CCD is also characterized by a root mean square (rms) readout noise, dark 
current, and quantum efficiency. 

The program has the option to change some of the parameters in the calcula- 
tions without restarting and completely re-typing all of the input data. 


USER INSTRUCTIONS 


The user is first asked to input the visual magnitude of the extended source 
and then to choose whether or not the spectral distribution of the source is uniform 
or blackbody. The cosmic background will also have the same type spectral distribu- 
tion. If a blackbody distribution is chosen, the user is asked to input the effective 
temperature (°K) of the source. The effective temperature for the cosmic back- 
ground will be requested later. 

Next the user is asked to input the lower and upper wavelengths (cm) and the 
number of wavelengths to consider. For Space Telescope, the lower and upper wave- 

-5 _it 

lengths are 10 cm and 11 x 10 cm. The maximum number of wavelengths that 
the program can handle is 21, which makes the wavelength increment 0.5 x 10 cm. 



For this case, the filter functions (0 < transmission < 1) and quantum efficiencies 
(electrons /photon) must be specified and read in (all 21 of them) in increments of 
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0.5 x 10 cm. The user is asked if there is a table of quantum efficiencies and 
filter functions. If the answer is yes, then two files of data must have been pre- 
viously prepared, which are now read into the program. One file contains the quan- 
tum efficiencies and the other contains the filter transmissions. Because the program 
integrates over wavelength, the first and last non-zero filter transmissions should be 
divided by two in order to use the trapezoidal rule for numerical integration. If the 
answer is no, the user is assumed to desire a system with 100 percent quantum 
efficiency and no filters. 

To obtain a monochromatic source and background, the lower wavelength will 
be chosen as the desired wavelength, if the number of wavelength calculations is set 
to one. The user is then asked to input the quantum efficiency and filter trans- 
mission for that wavelength. 

Next, the user is asked to input the diameter (cm) of the telescope entrance 
pupil and the focal ratio of the optical system. Then the visual magnitude of the 
cosmic background is requested. If a blackbody distribution was requested, an 
effective background temperature (°Kelvin) will also be requested. 

The user is then asked to specify the characteristics of the CCD pel array and 
the location of the extended source image within the array. To visualize the effects 
of the following inputs, see Figure 1 and consider the definitions listed below: 


XSW 

YSW 

XPW 

XPW 

XPD 

YPD 

XOF 

YOF 


X direction width of the extended source in arc sec 

Y direction width of the extended source in arc sec 
X direction width of the CCD pel in arc sec 

Y direction width of the CCD pel in arc sec 

X distance between pel centers in arc sec, YPD > XPW 

Y distance between pel centers in arc sec, YPD > YPW 

X distance offset in arc sec of extended source image from the origin, 
0 < XOF < XSW 

Y distance offset in arc sec of extended source image from the origin, 
0 < YOF < YSW . 


Dead space between the pels is created whenever XPD > XPW and/or YPD > YPW, i.e., 
the pel center separation is greater than the pel dimension in the same direction. 

The program computes the size of the array needed based upon the size and location 
of the extended source image. All parts of the source image fall either upon a dead 
space or a pel. 


After entering the x,y pel widths, x,y source image widths, x,y pel center 
separations, and the x,y source image offsets, the program computes a combined 
signal- to- noise ratio for an area that contains all pels effected by the extended source 
image. In the case of Figure 1, a combined signal- to-noise ratio would be computed 
for a three by three array. The user can obtain a S/N for a single pel completely 



filled by the portion of the source image by entering a "l," when the prompt "Enter 
Fraction of a Pel for S/N Per Pel" is received. To obtain the S/N for a pel partially 
filled by the portion of the extended source image, the user enters the fraction of 
the pel that is filled by the portion of the extended source. To complete the detec- 
tor characteristics, the user inputs the rms readout noise in the electrons per pel 
and the detector temperature in °Kelvin. For the particular CCD array used with 
the Space Telescope and the Wide Field /Planetary Camera, the readout noise has been 
suggested [1] to range from 13.9 to 17.8 electrons per pel rms. For the above CCD 
array, the maximum operating temperature [1] is expected to be around 178°K 
(-95°C), and the dark current in electrons per pel per second is computed from this 
input temperature. The user can input a mean dark current directly by first enter- 
ing a temperature less than 4°K. A prompt will then appear requesting the desired 
dark current. 

The next set of inputs is concerned with the output data that is desired. The 
user first specifies the observation start and end times in seconds and the number of 
time calculations. For example, to compute the S/N in hour intervals from 1 to 14 hr, 
the user would input a start time of 3600 sec, and an end time of 50400 sec and 14 
as the number of time calcualtions . To reverse the situation and determine the 
observation times needed to achieve a desired range of S/N, the user inputs a start 
and end S/N and the desired number of S/N calculations . The program will proceed 
to output the desired data. 

After the output is complete, the user may repeat the calculations for a source 
with a different visual magnitude and/or a background with a different visual magni- 
tude, without having to re-input all the previous input data. An additional repeat 
calculation that can be accomplished with or without the magnitude changes, is to 
change the characteristics of the detector. This option, however, requires that the 
user re-input all data starting with the detector characteristics and continuing through 
the number of S/N calculations. The program terminates when the user answers no 
to all the repeat calculation options. 


EQUATIONS USED 


The system of units used is the cgs system, except for the image source and 
pel dimensions, which are in arc-sec. The total energy emitted in the eye responsive 
spectral region per square arc-sec ( ® ) from an extended, astronomical source of 
visual magnitude m and received per square meter per second outside the Earth’s 
atmosphere is given by [2] 


I m * = 2.54 x 10 6 x 10 °* 4m lux/ GS (1) 

Equation (1) is converted from photometric to radiometric units by dividing by 
7.6x10 ^cm 

0.68 f K( A) dA lux per (ergs cm -2 sec -1 A A -1 ) , (2) 

3.8xl0 _5 cm 
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where A A is the wavelength interval (expressed in cm) and K( A) is the photopic eye 
response [2] given in Table 1. Dividing equation (1) by equation (2) is tantamount 
to assuming a uniform spectral distribution of the extended source and gives 

3.7353 x l(f 6 x l(f 0 - 4 m ergs cm -2 ® _1 sec' 1 A A' 1 

I m (A) = rg — (3) 

7.6x10 cm 

j K( A) dA 

3.8x10 ®cm 


If B( A ,T) is the blackbody distribution having an effective temperature T(°Kelvin), 
then equation (3) can be converted to a blackbody distribution in the following 
manner [ 3] : 




3.7353 x 10' 6 x 10 0,4 m B(A,T) ergs cm" 2 BO _1 sec" 1 A A -1 

7.6xl0~®cm 

J K(A) B(A,T) dA 
3.8x10 ®cm 


(4) 


Any type distributions can be used for B(A,T) in equation (4) since multiplying 

equation (4) by 0.68 K(A) and integrating from 3.8 x 10® cm to 7.6 x 10 ® cm 
reproduces equation (1). If B(A,t) were a uniform distribution, then equation (3) 
would be obtained again. Assume that a uniform distribution has been selected and 
that the following quantities have been specified: 


dft = area of the extended source in 0 

o 

A = area of the telescope entrance pupil (cm ) 
he /A = ergs /photon 

Q(A) = quantum efficiency in electrons per photon 
F(A) = optical filter transmission 


The quantum efficiencies [1] used in the output example are shown in Table 2. They 
include the combined effects of the optical telescope assembly and the wide field and 
planetary camera. The filter function [4] used in the output example is listed in 
Table 3. 

The rate at which electrons are liberated from the detector per wavelength 
interval at a chosen wavelength is given by 

S s (A) = - AI m (A)Q( A)F(A) electrons sec 1 A A 1 (5) 
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Equation (5) is called the monochromatic source signal current. Integrating equation 
(5) over the wavelength region influenced by the quantum efficiency and optical filter 
produces the polychromatic signal current. If there is a dead space between the 
detector pels, dfi is reduced to account for the dead space and any image offsets 
input by the user. The polychromatic source signal current is given by 


_ A(dn) 

s he 


X 2 

J I m ( A)Q( A)F( A) AdA electrons sec 


( 6 ) 


By the same analogy, the polychromatic, cosmic background, signal current is given 
by 


S 


b 


A(N p dfi') 

he 



AI m '(A)Q( A)F(A)dA electrons sec 


(7) 


where m' could be different visual magnitude than that of the source, dq’ is the area 
of single pel in ® , and N p is the number of pels that have been totally or partially 

illuminated by the extended source image. Thus, (N dfi’) is always equal to or 
greater than (dfi). p 

To compute a signal- to-noise ratio, the contributions to the signal term and the 
noise term must be defined. For the polychromatic case, for example, the signal is 
equation (6) multiplied by the time in seconds. If the incoming photon flux is 
assumed to be Poisson distributed, then the variance of the signal is also equal to 
the mean signal. If the signal to noise ratio is defined as the ratio of the mean to 
the standard deviation of the mean, then the source signal will contribute to the 
noise. Other contributions to the noise include the cosmic background signal, the 
rms readout noise in electrons per pel and the mean dark current in electrons per 
pel second. If all of the noise contributors are assumed to be statistically independ- 
ent, then the signal- to-noise ratio is given by 


(S/N) = 




(S g .t + S b .t+N p .R 2 +N p .D.t) 1/2 


( 8 ) 


where R is the rms readout noise, N is the total number of pels that are totally or 

ir 

partially illuminated by the extended source, D is the dark current and t is the time 
in seconds. As previously mentioned, the user has the option of inputing the dark 
current directly or inputing the detector temperature and having the dark current 
computed for the CCD array on Space Telescope for the Wide Field /Planetary Camera. 
The equations [5] for the dark current are given by 
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exp (-5802. 1Z) , 


(9) 


D = 29.3 x 10 9 x T 3/2 


where 


^ _ 1.1557 7. 021Txl0 * 

T 1108+ T ’ uu; 

where T is the detector temperature in °Kelvin. 

Equation (8) is a combined signal- to-noise ratio for all pels totally or partially 
illuminated by the extended source. To compute the signal to noise ratio for one pel 
that is totally or partially illuminated by the source, equation (8) is rewritten as 


(S/N)’ = 


(S s ’.t.f+S b .t+R 2 +D.t) 1/2 


(ID 


where f is the fraction of a pel illuminated by the source, and S ' is equation (6) 

s 

except that dn has been repalced with dfi’ as in equation (7). Equations (8) and 
(11) give the signal-to-noise ratios as a function of observation time. However, the 
time can be solved for in terms of the signal to noise ratio. Solving equation (8) 
for time gives 


2 <WV D) 

t = (s/Nr , - p - + 

2S S 


(S/N) 


4 <WV D > 2 2 

— / + (s/n) 2 n (r/s y 

4S 4 P s 


1/2 


sec 


( 12 ) 


Equation (11) can be similarly written by setting N =1 and replacing S with S’ f . 

p s s 
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TABLE 1. PHOTOPIC EYE RESPONSE 


(10 5 cm) 

K( A) 

A (10 5 cm) 

K( A) 

>-> 

/■— \ 

H- 1 

0 

l 

0 01 

3 

K( A) 

3.8 

0.00004 

5.1 

0.503 

6.4 

0.175 

3.9 

0.00012 

5.2 

0.71 

6.5 

0.107 

4.0 

0.0004 

5.3 

0.862 

6.6 

0.061 

4.1 

0.0012 

5.4 

0.954 

6.7 

0.032 

4.2 

0.004 

5.5 

0.995 

6.8 

0.017 

4.3 

0.0116 

5.6 

0.995 

6.9 

0.0082 

4.4 

0.023 

5.7 

0.952 

7.0 

0.0041 

4.5 

0.038 

5.8 

0.87 

7.1 

0.0021 

4.6 

0.06 

5.9 

0.757 

7.2 

0.00105 

4.7 

0.09 

6.0 

0.631 

7.3 

0.00051 

4.8 

0.139 

6.1 

0.503 

7.4 

0.00025 

4.9 

0.208 

6.2 

0.381 

7.5 

0.00012 

5.0 

0.323 

6.3 

0.265 

7.6 

0.00006 



TABLE 

2. QUANTUM EFFICIENCIES 


10 5 cm) 

Q(A) 

A (10 5 cm) 

Q( A) 

A (10 5 cm) 

Q( A) 

1.0 

0.0* 

4.5 

0.08* 

8.0 

0.12 

1.5 

0.028* 

5.0 

0.12 

8.5 

0.105* 

2.0 

0.05 

5.5 

0.15* 

9.0 

0.09 

2.5 

0.05 

6.0 

0.18 

9.5 

0.08* 

3.0 

0.05 

6.5 

0.18* 

10.0 

0.07 

3.5 

0.06 

7.0 

0.18 

10.5 

0.045* 

4.0 

0.04 

7.5 

0.15* 

11.0 

0.02 


*Interpolated Values 


TABLE 3. FILTER FUNCTION 


.0 5 cm) 

F(A) 

5.0 

0.278* 

5.5 

0.899 

6.0 

0.578 

6.5 

0.0835* 


♦The values were divided by 2 for using trapezoidal 
rule in numerical integration over wavelength. 
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SOURCE AND BACKGROUND 
ILLUMINATION 



BACKGROUND ILLUMINATION 
ONLY 



DEAD SPACE AND BACKGROUND 
ILLUMINATION 



DEAD SPACE, SOURCE AND 
BACKGROUND ILLUMINATION 


FIGURE 1. RELATIONSHIP OF EXTENDED SOURCE 
IMAGE TO CCD CONFIGURATION 
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APPENDIX A 


OUTPUT EXAMPLE 


ENTER EXTENDED SOURCE UISUAL MAGNITUDE 

72S . 

ENTEP CHOICE OF SPECTRAL DISTRIBUTION FOR SOURCE 
1-ELACKBODV OR O-FLAT 

70. 

ENTER LOUER AND UPPER UAUELENGTHS IN CM 

71. E-S.11.E-S 

ENTER NUMBER OF UAUELENGTH CALCULATIONS 


IS THERE A TABLE OF QUANTUM EFFICIENCIES?! YES* li NO-0) 
?1 . 

? IS THERE A TABLE OF FILTER FUNCTIONS?(YES-liNO-0) 

ENTER TELESCOPE ENTRANCE PUPIL DIAMETER! CM) AND 
SYSTEM FOCAL RATIO 

7240.. 30. 

ENTER COSMIC BACKGROUND UISUAL MAGNITUDE 
?S3. 

ENTER X,V WIDTHS (ARC SECONDS) OF SOURCE 
•’0.25,0.25 

ENTER X,V UIDTHSt ARC SECONDS) OF A PEL 

70.043.0. 043 

ENTER X, Y PEL CENTER SEPARATIONS! ARC SECONDS) 

70.043.0. 043 

ENTER X.Y DETECTOR OFFSET! ARC SECONDS) 

X,Y UALUES MUST BE -,>0,BUT < PEL UIDTHS 

70.004.0. 004 

ENTER FRACTION OF A PEL FILLED FOR !S'N) PER PEL 


ENTER RMS READOUT NOIS£!£LECTRONS/PEL > 

718. 

ENTER DETECTOR TEMPERATURE !KELUIN) 

7178. 

ENTER OBSERUATION START AND END TIMES ! SECONDS ) 

73600 . , 50400 * 

ENTER NUMBER! >1 ) OF TIME CALCULATIONS 
•’14. 

ENTER START AND END SIGNAL TO NOISE RATIOSiS/N) 

'’l. ,10. 

ENTER NUMBER! >1 ) OF (S/N) CALCULATIONS 

•’ 10 . 

CLEAR SCREEN AND RETURN FOR OUTPUT 
*? 
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79 . 


DO VOU UISH TO CHANCE BACKGROUND NACNITUDE?( YES-1, NO-O ) 
T0* 

DO VOU UISH TO CHANCE DETECTOR EFFECTS? (YES- 1; NO-0 ) 

?e. 

DID VOU HAKE ANY CHANGES? < YES- 1; NO-0 > 

?0 • 

ESTOP* 0 



oooooooooooooooo 


APPENDIX B 


PROGRAM LISTING 

07:29 AUG 10, '84 DC/EXTSRC. JAYRQE 

***** THIS PROGRAM COMPUTES THE SIGNAL TO NOISE RATIO ***** 
***** OUTSIDE THE EARTH'S ATMOSPHERE FOR A RECTANGULAR **** 
***** UNIFORM EXTENDED SOURCE OF VISUAL MAGNITUDE XSM ***** 
***** EMBEDDED IN A BACKGROUND OF VISUAL MAGNITUDE XMB **** 
***** AND FOR A RANGE OF OBSERVATION TIMES, THE CHAR- ***** 
***** ACTERIST ICS OF THE RECTANGULAR DETECTOR ARRAY ***** 
***** THAT ARE USED INCLUDE QUANTUM EFFICIENCY, PEL ***** 
***** DIMENSION, DEAD SPACE, RMS READOUT NOISE AND MEAN ***** 
***** DARK CURRENT. THE SOURCE AND BACKGROUND CAN HAVE ***** 
***** A FLAT SPECTRAL DISTRIBUTION OR A 8LACKB0DY DIS-***** 
***** TRIBUTION, IN WHICH CASE THE EFFECTIVE TEMPERA- ***** 
***** TURE OF THE SOURCE AND BACKGROUND MUST BE INPUT.***** 
***** THE OBSERVATION TELESCOPE/DETECTOR SYSTEM HAS AN***** 
***** ENTRANCE PUPIL DIAMETER DP AND A FOCAL RATIO FR,***** 
***** OPTICAL FILTERS CAN ALSO BE USED WITH FILTER ***** 
***** FUNCTIONS FF. ***** 

DIMENSION XK(40),QUEF(50),FF(50) 

WRITE ( 102,6000 ) 

6000 FORMATC' ENTER EXTENDED SOURCE VISUAL MAGNITUDE') 
READ(101,5000)XSM 
5000 FORMAT (20G) 

C ***** UNITS OF XSM ARE ( 1 0** ( -5 ) ) ERGS/CM*CM/SEC ***** 

XIMS=0, 37353* 1Q.0**(-0.4*XSM) 

WRITE(102,6004) 

6004 FORMAT!/,' ENTER CHOICE OF SPECTRAL DISTRIBUTION FOR SOURCE') 
WRITE(1Q2,6005) 

6005 FORMAT!' 1=BLACKBQDY OR OsFLAT ' ) 

DL A= 1 »E-6 

HC=1,9865E-16 
HCQK=1 .4388 
READ!lOt,5O00)KRESP 
IF (KRESP.EQ. 1 ) GO TO 100 

C ***** PHOTOPIC INTEGRAL OF K!L)DL FROM (3,8-7,6)E-5 ***** 

C ***** IS 1.0682E-5 CM. ***** 

XIMLS=XIMS/!1.0682*HC) 

GO TO 130 
100 SUMsO.O 

DO 110 IK=1,39 
110 READ!20,5000)XK!IK) 

WRITE! 102,6009) 

6009 FORMAT!/,' ENTER SOURCE EFFECTIVE TEMPERATURE !KELVIN) ' ) 

READ! 101, 5000) $E FT 

C ***** COMPUTE PHOTOPIC/BLACKBODY INTEGRAL (SIMPSON'S RULE) ***** 
X=3.8E-5 
Y=HCOK/(X*SEFT) 

SUMs$UM+XK(l)/((EXP(Y)-!,)*X**5) 

X=X+DLA 

DO 120 IK=2,38 
ICHEKe2*(IK/2) 

FAC *2, 
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IF!IK.EQ.ICHEK)FAC=4. 

Y=HCOK/!X*SEFT) 

SUM»SUM+FAC*XK(IK)/( (EXP(Y)-l , )*X**5) 
x=xVdla 
120 CONTINUE 

Y=HCOK/(X*SEFT) 

SIJMT=SUM + XK!39)/! ! EXP ! Y ) - 1 , ) *X**5 ) 

XIMLS=XIMS*SUMT*DLA/!3,*HC) 

130 WRITE! 102,6013) 

6013 FORMAT!/,' ENTER LOWER AND UPPER WAVELENGTHS IN CM') 

READ! 101,5000) WAVL,WAVH 

WAV=WAVL 

WRITE (102,6017) 

6017 FORMAT!/,' ENTER NUMBER OF WAVELENGTH CALCULATIONS') 
READ!101,5000)NINC 
IF!NINC.EQ,l)GO TO 180 
XsNINC-1 

DL=!WAVH-WAVL)/X 
DO 140 1=1,21 
QUEF ( I ) = 1 , 

FF ( I ) = i , 

140 CONTINUE 
FF ! 1 )=0,5 
FF!21)=0.5 
WRITE! 102,6021) 

6021 FORMAT!/,' IS THERE A TABLE OF QUANTUM EFF IC IENC IES? ! YES= 1 1 N0 = 0 ) ' ) 
READI101, 5000)1 RE SP 
IF!IRESP.EQ.O)GO TO 160 
DO 150 1=1 , NINC 
150 READ!22,5000)QUEF!I) 

160 WRITE! 102,6025) 

6025 FORMAT!/,' IS THERE A TABLE OF FILTER FUNCT IONS? ! YES= 1 * NQ = 0 ) ' ) 
READ!10i,5000)IRESP 
IF ! IRESP.EQ, 0)G0 TO 200 
DO 170 1=1, NINC 
170 READ(24,5000)FF!I) 

GO TO 200 

180 WRITE ! 102,6029) 

6029 FORMAT!/,' ENTER MONOCHROMATIC QUANTUM EFFICIENCY') 
READ!101,5000)9UEF!1) 

WRITE!102,6033) 

6033 FORMAT!/, 'ENTER MONOCHROMATIC FILTER FUNCTION') 

READ!10l,5000)FF!l) 

IF (KRESP»EQ* 1 )GO TO 190 
SUM=XIMLS*WAVL*QUEF!1 )*FF!1 ) 

GO TO 260 

190 Y=HCQK/(WAVL*SEFT) 

SUM=XIMLS*QUEF (1)*FF!1)/!(EXP(Y)«1, )*WAVL**4) 

GO TO 260 
200 MINC*NINC-1 

IFIKRESP.EQ.l )G0 TO 230 
SUM=WAV*QUEF!1)*FF!1) 


13 



IF ( M I NC ,EQ.2)G0 TO 220 
DO 210 I-2/MINC 
WAV=WAV+PL 
PROD=QUEFCI)*FFCI) 

IF(PR0D.LT.1,E-4)G0 TO 210 

sum=sum+wav*prod 

210 CONTINUE 
220 WAV=WAV+DL 

SUM=SUM+WAV*QUEF(NINC)*FF(NINC) 

SUM=SUM*DL*XIMLS 
GO TO 260 

230 Y=HCOK/ ( WAV*SEFT) 

SUM=QUEF(l)*FF(l)/( (EXP(Y)-l . )*WAV**4) 

IF(MINC.EQ.2)G0 TO 250 
DO 2^40 1=2 # MINC 
WAVsWAVtOL 
PROD=QUEF (I)*FF(I) 

IFCPROD.LT. l.E-4)G0 TO 240 
Y=HCOK/(WAV*SEFT) 

SUM=SUM+PR0D/((EXP(Y)-1.)*WAV**4) 

240 CONTINUE 
250 WAVsWAVfOL 

Y=HCOK/(WAV*SEFT) 

SUM=SUM+QUEF (NINC) *FF (NINC)/((EXP(Y)*1.)*WAV**4) 
SUM=SUM*DL*XIMLS 
260 WRITE(102,6037) 

6037 FORMAT (/# 'ENTER TELESCOPE ENTRANCE PUPIL DIAMETER(CM) AND') 
WRITE(1Q2#6Q38) 

6038 FORMAT ( 'SYSTEM FOCAL RATIO') 

READ(101,50Q0)DP,FR 
SIGNL=0.7854*DP*DP*SUM 
WRITE ( 102#6041 ) 

6041 FORMAT!/, 'ENTER COSMIC BACKGROUND VISUAL MAGNITUDE') 
READ(101,5000)XM8 
XIM8=0. 37353* 10.0**(-0,4*XMB-5, ) 

WAV=WAVL 

IF (KRESP.EG. 1 )GO TO 270 
XIMLB=(XIMB*10,0**5)/C1,Q682*HC) 

GO TO 285 

270 WRITE (102, 6049) 

6049 FORMAT (/# 'ENTER BACKGROUND EFFECTIVE TEMPERATURE (KELVIN) ' ) 
READC 101 #5000) BEFT 
X=3.8E-5 
Y=HCQK/(X*BEFT) 

SUM=XKCl)/C(EXP(Y)-l,)*X**5) 

XSX.+DLA 
DO 280 IK=2# 38 
ICHEK=2* ( IK/2) 

FAC=2, 

IFCIK.EQ* ICHEK)FAC=4. 

Y=HCOK/(X*BEFT) 

SUM=SUM+FAC*XK CIK)/ ( (EXP(Y.)-l • )*X**5) 
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XSX+DLA 
200 CONTINUE 

YaHCOK/(X*BEFT) 

SUM=SUM+XK(39)/((EXP(Y)-1.)*X**5) 

XIML8aXIMB*3./(SUM*DtA*HC) 

285 IF (NINC.GT# 1 J GO TO 300 
IF (KRESP.EQ. 1 )G0 TO 290 
SUM=XIMLB*WAV*QUEF(l)*FF(l) 

GO TO 360 

290 Y=HCOK/(WAV*BEFT) 

SUM=XIML8*QUEF(1)*FF(1)/((EXP(Y)-1.)*WAV**9) 

GO TO 360 

300 IF (KRESP,£Q. 1 )G0 TO 330 
$UM=WAV*QUEF ( 1 )*FF ( 1 ) 

IF (MINC.EQ,2)G0 TO 320 
DO 310 I=2,MINC 
WAVsWAV+DL 
PROD=QUEF(I)*FFCI) 

IF(PR0D.LT.1.E-9)G0 TO 310 
SUMaSUM+WAV*PROD 
310 CONTINUE 
320 WAV=WAV+DL 

SUM=SUM+WAV*QUEF(NINC)*FF(NINC) 

SUM»SUM*DL*XIMLB 
GO TO 360 

330 Y=HCOK/(WAV*BEFT) 

SUMsQUEF(l)*FF(l)/((EXP(Y)-l.)*WAV**9) 

IF (WINC.EQ.2JG0 TO 350 
DO 390 I=2,MINC 
W A V=WA V+DL 
PROD*QUEF(I)*FF(I) 

IF(PR0D.LT.1.E-9)GQ TO 390 
Y=HCOK/(WAV*BEFT) 

SUM=SUM+PR0D/((EXP(Y)-1,)*WAV**9) 

390 CONTINUE 
350 WAV=WAV+DL 

Y=HCOK/(WAV*BEFT) 

8UM=SUM+QUEF(NINC)*FF (NINC)/ ( (EXP (Y)«l . )*WAV**9) 
SUM=SUM*OL*XIMLB 

360 BGRND=0.7859*DP*DP*SUM 

361 WRITE C 102*6053) 

6053 FORMATC/# 'ENTER X,Y WIDTHS(ARC SECONDS) OF SOURCE') 
READ(101,5000)XSW,YSW 
WRITE (102,6057) 

6057 FORMAT (/* 'ENTER X,Y WIDTHS(ARC SECONDS) OF A PEL') 
READ(101,500Q)XPW,YPW 
WRI TE ( 1 02/606 1 ) 

6061 FQRMATC/r'ENTER X,Y PEL CENTER SEPARATIQNStARC SECONDS)') 
R£AD(101,50QO)XPO,YPP 
WR ITE ( 1 02,6065 ) 

6065 FORMATC/, 'ENTER X,Y DETECTOR OFFSET (ARC SECONDS)') 

WRI TE ( 1 02,6069) 



6069 FORMAT ! 'X, Y VALUES MUST BE =,>0,BUT < PEL WIDTHS') 
READ(101,5000)X0F # YQF 
XN=!XOF+XSW)/XPD 
NX = XN 
MX = NX 
XNX=NX 

IF(XN,GT.XNX)MX=NX+1 
XMsMX 
p R J X s 0 

IFCXPdIlT.XPWJGQ TO 365 

PRTXs(XQF+XSW-XNX*XPD-XPW)/(XPD-XPW) 

IF(PRTX.LE.0)PRTX=0. 

365 YN=(YQF+YSW)/YPD 
NY = YN 
MY = NY 
YNYsNY 

IF(YN,GT.YNY)MY=NY+1 

YM=MY 

PRTY=0. 

IF ! YPD.LE. YPW)GO TO 367 
PRTY*(YOF+YSW-YNY*YPD-YPW)/(YPD*YPW) 

IF(PRTY.LE.Q.)PRTYsQ. 

367 A$=XSW*YSW+!XPD-XPW)*(YPD-YPW)*(XNX*YNY+PRTX*YNY+PRTY*XNX 
&+PRTX*PRTY) 

AS=AS-YSW*!XPD*XPW)*!XNX+PRTX)-XSW*!YPD-YPW)*(YNY+PRTY) 

AB=XPW*YPW*XM*YM 

TOTS*SIGNL*AS 

TQTB»BGRND*AB 

WRITE!102,6073) 

6073 FORMAT !/, 'ENTER FRACTION OF A PEL FILLED FOR (S/N) PER PEL') 
READ(1O1,5O0O)FRAC 
PELS=SIGNL*FRAC*XPW*YPW 
PEL8*BGRND*XPW*YPW 
WRITE! 102,6077 ) 

6077 FORMAT!/, 'ENTER RMS READOUT NOI SE (ELECTRONS/PEL )' ) 
READ(101,5QQ0)R0N 
RN=RON*RON 
WRITE!102,6081) 

6081 FORMAT!/, 'ENTER DETECTOR TEMPERATUREIKELVIN) ' ) 

READ ! 101,50 00 )DTEM 
IF!DTEM.LT.O.)GO TO 370 

FACs-6705.A87/DTEM+A,0737*DTEM/(n08,+DTEM) 

FAC=EXP!FAC)*!29,3E+9) 

DC=FAC*DTEM**l ,5 
GO TO 380 

370 WRITE! 102, 6085) 

6085 FORMAT!/, 'ENTER MEAN DARK CURRENT ! ELECTRONS/ PEL /SECOND) ' ) 
READU01,S000)DC 
380 WRITE ! 102,6089) 

6089 FORMAT!/, 'ENTER OBSERVATION START AND END TIMES(SECONDS) ') 
READ!101,5000)STRT,ENDT 
WRITE !102, 6093) 
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6093 FORMAT (/, 'ENTER NUM8EROD OF TIME CALCULATIONS') 

READ(101,5000)NT 

X=NT-1 

DT=(ENDT-STRT)/X 

WRITE(102,6094) 

6094 FORMAT (/, 'ENTER START AND END SIGNAL TO NOISE RATI0S(S/N) ' ) 
READ(101,5000)SNST,SNND 

WRITE ( 102,6095) 

6095 FORMAT!/, 'ENTER NUMBER C > 1 ) OF (S/N) CALCULATIONS') 

READ(101,5000)NS 

V=NS-1 

DNs(SNND-SNST)/Y 

WRITE(lO2,6096> 

6096 FORMAT!/, 'CLEAR SCREEN AND RETURN FOR OUTPUT') 

READdOl ,5000)DUMMY 

445 WRITE(102,6Q97)DP 

6097 FORMAT!/, 'TELESCOPE ENTRANCE PUPIL DIAMETER ',F7.3,' CM') 
WRITE(102,6101)FR 

6101 FORMAT! 'OPTICAL SYSTEM FOCAL RATIO ',F5.2) 

WRITE(102,6105) 

6105 FORMAT!/, 'SOURCE CHARACTERISTICS') 

IF (KRESP.NE. 1 )GG TO 390 
WRITE!102,6109)SIGNL 

6109 FORMAT! IX, 'PLANC DISTRIBUTION', E15, 8, ' DETECT. ELEC/SEC/ARCSEC**2 ' ) 
WRITE!102,6113)SEFT 

6113 F0RMATI5X, 'EFFECTIVE TEMPERATURE ',F9.3,' DEGREES KELVIN') 

GO TO 400 

390 WRITE!102,6117)SIGNL 

6117 FORMATI1X, 'FLAT DISTRIBUTION ', El 5. 8 , ' 0ETECT.ELEC/SEC/ARCSEC**2') 
400 WRITE!102,6121)XSM 

6121 F0RMATI5X, 'EXTENDED SOURCE MAGNITUDE ',F5.2) 

WRITE! 102, 61 22 )XI MS 

6122 FORMAT(5X,'ERGS/SEC/CM**2/ARCSEC**2 OUTSIDE ATMOSPHERE ',E15.P) 
WRITE!102,6125)XSW,YSW 

6125 F0RMAT!5X, 'SOURCE DIMENSIONS! ARC SECONDS ): X= ', F0 . 4, ' Y=',F8.4) 
WRITE ! 1 02, 6129) 

6129 FORMAT!/, 'BACKBGROtJND CHARACTERISTICS') 

IF (KRESP.NE. 1 )G0 TO 410 
WRITE(l02,bl09)BGRND 
WRITE!102,6U3)BEFT 
GO TO 420 

410 WRITE{102,6tl7)BGRND 
420 WRITE!102,6145)XMB 

6145 FORMATI5X, 'BACKGROUND MAGNITUDE ',F5.P) 

WRITE!102,6122)XIMB 
WRITE (102,6149) 

6149 FORMAT!/, 'SYSTEM SPECTRAL CHARACTERISTICS') 

WRITE! 102,6153) 

6153 FORMAT (4X, 'WAVELENGTH (CM)', 5X, 'QUAN.EFFIC. ',3X, 'FILTER FUNCTION') 
WAV=WAVL“DL 
DO 430 1=1 , NINC 
WAV=WAV+DL 
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PRQD=QUEF(I)*FF(I) 

IF(PR0D.LT.1.E-4)GQ TO 430 
WRITE(102,6157)WAV,QUEF(I),FF(I) 

6157 FQRMAT(3C2X,E15.0)) 

430 CONTINUE 

WRITE ( 1 02/6161) 

6161 FORMAT C/ r 'DETECTOR CHARACTERISTICS') 

WRITE C 1 02,6165) XPW,YPW 

6165 FORMAT (5X, 'PEL DIMENSIONSCARC SECONDS) s X= ', F8. a , ' Y=',F8.4) 
WRITE(102,6169)XPD,YPD 

6169 FORMAT (5X, 'PEL CENTER SEPARATIONCARC SECONDS ): X= ', F8. 4 , ' Y=',F8.4) 
WRITE(102,6173)XOF,YOF 

6173 F0RMAT(5X, 'DETECTOR QFFSETCARC SECONDS) l X= ', F8, 4, ' Y=',F8.4) 
WRITE(102,6177)MX,MY 

6177 F0RMATC5X, 'ARRAY SIZE EFFECTED BY SOURCE : Xs ', 12 , ' BY Y=',I2) 

WRITE (102,6101 IRON 

6181 FORMAT (5X, 'RMS READOUT NOISE (ELECTRONS/PEL ):', F5. 2 ) 

IFCDTEM.LT. 4. )G0 TO 440 
WRITE Cl 02,61 85 )DTEM 

6185 F0RMATC5X, 'DETECTOR TEMPERATURE (KELVIN) :',F8.4) 

440 WRITE(102,6189)DC 

6189 FORMATCSX, 'DARK CURRENT (ELECTRONS/PEL/SECOND) : ',E15,8) 

WRITE C 1 02,6193) 

6193 FORMATC/, 'SIGNAL TO NOISE RATIO AS A FUNCTION OF OBSERVATION', 

&' TIME') 

WRITEC102/6197) 

6197 FORMAT ( 'PER ARRAY (S/N) AND PER FRACTION OF A PEL CS/N/P ), WHERE ' ) 
WRITEC102,6201)FRAC 

6201 FORMAT ( 'THE FRACTION IS ',F7,4) 

WRI TE C 1 02,6205) 

6205 FORMAT C4X, 'TIME (SECONDS) ',2X,' (S/N) ',2X,' (S/M/P)') 

PRON=XM*YM*RN 
PDC=XM*YM*DC 
T=STRT-DT 
DO 450 1=1, NT 
T=T+DT 
XNUM*TOTS*T 

XDNOMsXNUM+TOTB*T+PRON+PDC*T 

SNTsXNUM/SQRTCXDNOM) 

XNUMspELS*T 

XDNOMsXNUMtPELB*T+RN+DC*T 
SNP=XNUM/SQRT (XDNOM) 

WRITE(102,6157)T,SNT,SNP 
450 CONTINUE 

WRITE C 102,6206) 

620.6 FORMAT C/,5X, ' CS/N) ',4X, 'SECONDS (ARRAY) ',3X, 'SECONDS (PIXEL) ' ) 
SN=SNST-DN 
DO 455 1=1, NS 
SNsSN+DN 

A=SN*$N/(T0TS*TQTS) 

8=A*(TOTS+TOTB+PDC)*0.5 

C=A*PR0N 
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T1=B+SGRT(B*8+C) 

AsSN*SN/CPELS*PELS) 

B=A*!PElS+PELB+DC)*0,5 

C=A*RN 

T2=B+SQRT (8*8+C ) 

WRIT£U02,6207)SN,T1,T2 

6207 F0RMAT(5X,F5.2,2!2X,E15.8)) 

455 CONTINUE 

WRITE ( 1 02*6209) 

6209 FORMAT ! / , '00 YOU WISH TO CHANGE SOURCE MAGNI TUOE? ( YES = 1 j N0 = 0 ) ' ) 
READ! 101 ,5000) JRESP 
IF! JRESP.EQ.O)GO TO 460 
WR ITE ! 1 02 , 62 1 3 ) 

6213 FORMAT !/, 'ENTER NEW SOURCE MAGNITUDE') 

READ(101,5000)XSM 

XNEWsO, 37353* 10,Q**!-0,4*XSM) 

XNEWsXNEW/XIMS 
SIGNUsSlGNU*XNEW 
TOTS«TOTS*XNEW 
PELS=PELS*XNEW 
460 WRITE (1 02, 621 7) 

6217 FORMAT!/, 'DO YOU WISH TO CHANGE BACKGROUND MAGNI TUDE? ( YES = 1 ?' , 
&'NO=0) ') 

READC1O1,5QO0)JRESP 
IF!JRESP.EQ.Q)GQ TO 470 
WRITE(1Q2,6221 ) 

6221 FORMAT!/, 'ENTER NEW BACKGROUND MAGNITUDE') 

READ!101,5000)XMB 

XNEWs0.37353*10.0**!-0,4*XMB) 

XNEWsXNEW/XIMB 
BGRN0=8GRND*XNEW 
T0TB=T0T8*XNEW 
PEL8=PELB*XNEW 
470 WRITE! 102,6225) 

6225 FORMAT!/, 'DO YOU WISH TO CHANGE DETECTOR EFFECTS? ! YES= 1 ; N0=0 )' ) 
READOOl ,5000) JRESP 
IF (JRESP.EQ. 1 )G0 TO 361 
WRITE! 102,6229) 

6229 FORMAT!/, 'DID YOU MAKE ANY CHANGES? ( YES= 1 ; N0 = 0 )' ) 
READ!101,5000)JRESP 
IF (JRESP.EQ. l)GO TO 445 
STOP 
END 
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